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Abstract 

The study of the Dyson-Schwinger equations of Landau gauge Yang-Mills theory has revealed 
two types of solutions for the gluon and ghost propagators, with a scaling and a massive (decou- 
pling) behavior in the extreme infrared, respectively. We show that both types of solutions are 
quantitatively reproduced by applying renormalization group equations of Callan-Symanzik type in 
an epsilon expansion to the infrared limit of Landau gauge Yang-Mills theory when a mass term for 
the gluons is added to the action. Only the decoupling solution corresponds to an infrared-stable 
fixed point in three and four space-time dimensions and is hence expected to be physically realized, 
in agreement with the results of recent lattice calculations. 
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The infrared (IR) behavior of gluon and ghost propagators in SU{N) Yang-Mills theory 
in Landau gauge has been a subject of intense debate over the last decade. The first results 
that have been obtained in the extreme IR region, for momentum scales p <^ Aqcd, came 
from the solution of an approximation to the full Dyson-Schwinger (DS) equations in this 
theory l|]. They are nowadays known as scaling solutions, showing power-law behavior of 
the propagators in the IR. 

For a quantitative description of these results, we introduce some notation. The gluon 
and ghost propagators in D Euclidean space-time dimensions are parameterized as 



9 I "A''^ 9 
P \ P / 

{np)c\-q)) = ^ 5'^^(2vr)^5(p - q) . (1) 

The scaling solutions show a power behavior of the dressing functions in the IR, 

G{p') oc , oc ^ (2) 

The anomalous dimensions oip fulfill a sum rule, 

D-A 

aa + 2aF = — - — • (3) 

The approximate DS equations can be solved analytically in the IR regime . There 
are two consistent solutions, one of them giving 

«Fp) = ^, ao{D) = -^, (4) 

for 2 < D < 4. The IR exponents of the other solution are obtained, for 2 < < 4, in 
general numerically from a transcendental equation. However, a very good approximation 
(within 2% error) is given by 

«Fp) = ^, aaP) = -^^. (5) 

Note that, within the range of dimensions indicated, for both solutions ap > and ac < — 1 
hold, which means that the ghost propagator is enhanced in the IR with respect to the tree- 
level propagator, while the gluon propagator tends to zero for p — > 0. As a consequence, 
in the IR the diagrams with the greatest number of ghost propagators dominate over the 
others, a property called ghost dominance. 



Of particular interest is the IR behavior of a dimensionless running couphng constant 



defined from the ghost-gluon vertex. Taylor's well-known argument ^ implies that in the 
Landau gauge at the symmetrical point the proper ghost-gluon vertex receives only finite 
and scale- independent quantum corrections, to all perturbative orders. This fact motivates 
the definition 

gUp') = ip'r-'^^' G{p')F\p')g' (6) 

of the dimensionless renormalized coupling constant in terms of the bare coupling constant 
g. As a consequence of the sum rule ([3]), for the scaling solutions g^ tends toward a nonzero 
constant in the limit p — )■ 0. The value of gniO) can be obtained from the IR solutions of 
the DS equations. 

Several years after the discovery of the scaling solutions, however, another type of solu- 
;ions of the same approximate DS equations have been found, called decoupling solutions 
5|-[7|. The latter solutions are characterized hj ap = and ac = —1, for any dimension 
D in the range 2 < D < A. In particular, for these solutions the sum rule ([3]) does not 
hold. Qualitatively, in the decoupling solutions the ghost propagator is enhanced in the IR 
at most by a finite factor, while the gluon propagator tends towards a nonzero constant 
for p — )■ 0, similar to the propagator in a massive theory. Consequently, the dimensionless 
running coupling constant defined in Eq. ([6]) vanishes in the IR limit, contrary to the scaling 
solutions. 

Ever since the scaling solutions have been found, lattice calculations in the Landau gauge 
have focused on the extreme IR regime. The subject has been controversial over many 
years, but as for today there is quite convincing evidence that these numerical simulations 



find gluon and ghost propagators of the decoupling type in dimensions D = 3,4 [sHld] and 
scaling solutions in D = 2 

In the present paper, we show that a renormalization group (RG) analysis in an epsilon 
expansion is capable of quantitatively reproducing both the scaling and the decoupling so- 
lutions in the extreme IR, and furthermore of determining the IR-stability of the solutions. 
Only solutions corresponding to IR-stable fixed points are expected to be physically real- 
ized. Of course, it is well known that in the standard formulation of the theory the only 
perturbatively accesible fixed point is the trivial one. It is UV-attractive and IR-repulsive 
and can therefore not be used for a description of the IR physics. We here introduce a 
minimal modification of standard Yang-Mills theory by adding a mass term for the gluons 



as in Ref. 



12| , see also [3] • The gluon mass is generically expected to be of the order of the 
QCD scale Aqcd- It is well known that this modification does not affect the UV behavior 
of the theory. In particular, the beta function in this aymptotic region is unchanged and 
asymptotic freedom is maintained. A gluon mass term has been considered in many contexts 
before, often on phenomenological grounds. 

Our main motivation for the introduction of a gluon mass term, however, comes from new 
theoretical results: the quantization of Yang- Mills theory in the Landau gauge is complicated 
by the existence of gauge copies, and the corresponding restriction of the functional integral 
over the gauge fields to the first Gribov region breaks the BRST invariance of the effective 



action 
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15| . From the RG viewpoint, the appearance of a gluon mass term is a natural 



consequence of the breaking of BRST invariance. Indeed, it has been shown in Ref. [15| that 
a gluon mass term arises as a result of the formation of condensates of auxiliary fields that 
are introduced in order to implement the restriction to the Gribov region by adding local 
operators to the action. In the present analysis, we will not take any additional auxiliary 
fields into account. 

In the following, we will hence consider the D-dimensional Euclidean action 

S = Jd^x (^i + 1 mM» + d.c^Dfc' + zB'^d.A'l^ (7) 

with the gluon mass m, as the classical action for SU{N) Yang-Mills theory in the Landau 
gauge. The Nakanishi-Lautrup field is used to conveniently implement the restriction 
to transverse gauge fields, d^A'^ = 0. As we shall now argue, the existence of a gluon mass 
term profoundly changes the RG analysis of the theory in the IR. 

Naively, for momenta p with ^ m^, the mass term in Eq. ([7]) dominates the contri- 
butions to the action that are quadratic in the gluon field. In fact, this is also true for the 
loop corrections. As a relevant example, consider the one-loop contribution to the ghost 
self-energy in dimensions 2 < D < 4: it is straightforward to show that one may replace the 
tree-level gluon propagator corresponding to the action ([7]) inside the loop diagram with 

^{i..-'-f)i'WS(p-,) (8) 

and obtain exactly the same result in the IR limit, with corrections suppressed by powers 
of p^/m^, after renormalization with an appropriate normalization condition for the term 
cx in the ghost two-point function. For the description of the extreme IR regime, we will 



hence in the following reduce the complete term quadratic in the gluon field in the action 
([7]) to the mass term. 

We initiate the RG analysis by considering the limit g ^ of the approximate IR-form 
of the action which then reduces to 

5°, = j d^'x Q A" + d^c'^d^c^ + iB-d^AI^ . (9) 

Under a rescaling of the space-time coordinate x x/s with s > 1 [to be understood as one 
step in a Wilsonian RG transformation for the free field theory defined by the action (Q], 
the IR-action is invariant provided that the fields transform as 

(m^ being constant under the rescaling). Note that Eq. f ITU]) now defines the canonical or 
scaling (mass) dimensions of the fields which differ, in the case of the gluon field and the 
Nakanishi-Lautrup field, from the usual canonical dimensions which refer to the Gaussian 
fixed point. To reconcile the scaling dimension of with the usual canonical dimension, 
we may absorb one power of the mass in the gluon field to define 

A^^ix) = mA;{x) (11) 

and consider as the fundamental field. Similarly, we define B°-{x) = B"-{x)/m. As far as 
the gluon field is concerned, the fixed point action corresponds to the high-temperature 
fixed point rather than the usual Gaussian fixed point [loi]. It is natural to associate the 
high-temperature fixed point with a stochastic vacuum state. 

We will now reintroduce the coupling terms appearing in the Yang-Mills action. An 
important consequence of the scaling dimensions (fTOj) is that the ghost-gluon, three-gluon, 
and four-gluon couphngs all receive different scaling dimensions: 

gcAc s^-^^/^^ gsAc , ^ ^ _^ s-^gA4 , (12) 

so that all the coupling constants are perturbatively irrelevant [around the IR fixed point 
described by the action ([9])], except for gcAc at dimensions D < 2. In other words, D = 2 is 
the upper critical dimension of the theory, in agreement with an analysis of the IR behavior 
of higher orders in perturbation theory pJi] . Note that the scaling dimensions (|T2l) are simply 

5 



the mass dimensions of the couphng constants when one replaces the gluon field with the 
rescaled field according to Eq. ffTTl) . 

The theory we will hence be considering in the following is defined by the free action 
plus the sole interaction term 



(13) 



where we have redefined the coupling constant as ^ = gcAc/^- K will be interesting to 
consider this specific theory for dimensions D > 2, even though the only IR-stable fixed 
point for D > 2 will turn out to be the trivial one (Q, as expected. We will employ an 
epsilon expansion around D = 2, i.e., we will consider dimensions D = 2 + e. 

Note that neglecting the three- and four-gluon couplings (with the argument that these 
couplings are perturbatively irrelevant for any positive dimension) precisely corresponds 
to ghost dominance because then only the diagrams with the greatest number of ghost 
propagators remain. Ghost dominance can hence be considered as a consequence of the 
change in the scaling dimension of the gluon field which, in turn, is a consequence of the 
appearance of a gluon mass term. 

For the RG analysis of the interacting theory we calculate the one-loop self-energies, 
first for the ghost fields. From the action plus the interaction term fll3p one obtains 
(suppressing the color factor) 



1 Ng 



^2 



2 47r 



■P 



2 V 
^ ' + 7e - ln(47r) + In - 

K 



(14) 



at D = 2 + e, where we have introduced a dimensionless coupling constant g via 



~9' = {^'r^'f 



(15) 



[cf. Eq. ( 1T2|) ]. For the gluon self-energy, on the other hand, one obtains in D = 2 + e 
dimensions. 



1 Nf 



2 4n 



- + 7^-ln(47r)+ln^-2')5,. + 2^ 



(16) 
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Note that the gluon self-energy contains a longitudinal contribution which makes the break- 
ing of BRST invariance as a result of ghost dominance explicit. In the following we only 
consider the transverse part of the self-energy (since the gluon propagator stays transverse) 
which is given by the coefficient of 6^1, in Eq. f ITB]) . We also emphasize that the counterterm 
for the self-energy fll6p can be chosen local. 

We now renormalize the theory in the standard way, rescaling the fields with field renor- 
malization constants Za and Zc, and impose the following normalization conditions on the 
propagators of the renormalized fields: 

(4(p)4(-?)> , , = ^ {o\v)A-^)) , ^ = ^5^^(27r)^5(p-g). (17) 

Upon substituting the expressions (JH]), f lT6|) for the self-energies, the renormalization con- 
stants become functions of the renormalization scale /i. The resulting one-loop anomalous 
dimensions read 

The /^-independence of the bare propagators implies, in general, that 

/^'^ <4(p)4(-?)> = -7c (4(p)4(-?)> • (19) 

For the definition of the (dimensionless) renormalized coupling constant we consider the 
renormalized proper ghost-gluon vertex at the symmetric point, 

= /^-^/' 9Rif^)r'''Wp (2vr)^5(j9 + q + r). (20) 

(We have suppressed the Lorentz and color indices on the proper ghost-gluon vertex Fg^^ for 
readability.) Simple power counting shows that the loop corrections to Fg^^ at the symmetric 
point are /i-independent (Taylor's theorem). Consequently, one finds for the beta function 
to one-loop order, 

P(Sr)-^'-4,Sr-1;sJ-,-1^)^ (21) 



For e > (D > 2), we find a trivial IR-stable fixed point (as expected) and a nontrivial 
IR- unstable fixed point. Although the latter is irrelevant for the IR physics, for the sake 
of comparison with the results of the DS analysis it is interesting to determine the IR 
behavior of the propagators corresponding to this fixed point. Substituting the unstable 
fixed-point value of in Eq. fllSp and integrating the differential equations (IT^ for the 
renormalized propagators with these fixed point values for the anomalous dimensions and 
the normalization conditions f|T71) as initial conditions yields 



<4(p)4(-?)> = ^ [^y'n^^)''SiP-Q)- (22) 

The IR-limit of the propagators is hence of the form ([2]), and the exponents coincide exactly 
with the scaling solution @ of the DS equations found in Refs. {2, 3|. 

The other fixed point in Eq. (|2T1) . which is IR-stable and hence physically relevant for 
the description of the IR regime of the theory in D > 2 dimensions, is trivial in accord 
with our identification of D = 2 as the upper critical dimension. Apart from a vanishing 
renormalized coupling constant, one hence finds the "classical" behavior of the propagators 
in the IR limit, which means in particular that (the transverse part of) the gluon propagator 
tends toward a non- vanishing constant for zero momentum [cf. Eq. (|9])]. The trivial fixed 
point then corresponds to the decoupling solution js-T] of the DS equations described in the 
paragraph after Eq. ([6]). The triviality of the fixed point justifies the use of perturbation 
theory for the decoupling solution [13]. 

We can obtain a more precise description of the IR behavior of the propagators and the 
renormalized coupling constant by analyzing the approach to the (stable) IR fixed point as 
the momentum scale tends to zero. Namely, integrating the differential equation f l2T|) for 
gRifi) gives 

Air l + (/i7A2)^/2 ' ^ ^ 

where A is the (undetermined) scale with 

"halfway" between the two fixed points. Substituting the expression for gR^fi) from Eq. 
(!23|l for g in the anomalous dimensions ( 1181) and using the result in the integration of the 
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differential equations (fT9|l for the renormalized propagators [together with Eq. ( !T7|) ] leads to 

<4(p)4(-^)> = ^ it(pV%^-/^' ^'^^(2vr)-5(p - g) . (25) 

Also note that, from Eq. ([6]), gjiip^) = = P^) (except for a constant factor from the 

finite renormalization of the ghost-gluon vertex). 

"or the propagators is at least qualitatively confirmed by calculations on 
8 



The result ([25]) 
very large lattices 



10l |. In dimension D = 3, the linear rise of the gluon propagator with 
\p\ from a constant finite value in the IR is clearly seen jl7| (incidentally, this is very different 
from the usual behavior of a massive propagator). In Z) = 4 dimensions, a (approximately) 
quadratic rise of the gluon propagator with \p\ in the IR as predicted by Eq. ( 125|) is found 
in lattice simulations with lattice parameter /3 = 0, i.e., without any contribution from the 
original Yang- Mills action |18j. In order to compare our results in four dimensions to the 
case of finite /3, we have to reintroduce the term (A^p^A^) in the action, since the latter 
term is of the same order in (jp /m?) as the correction to the "classical" behavior generated 
by the application of the RG (at e = 2). Given that there are no loop corrections to the 
(y4jjp^y4^)-term at one-loop order, the resulting gluon propagator is precisely of the form 
obtained in the refined Gribov-Zwanziger scenario 15|. 

We also mention that a semi-quantitative reproduction of the lattice results (for the de- 
coupling solution) over the whole momentum range was achieved in D = 3 and 4 dimensions 
in one of the renormalization schemes proposed in Ref. 13|]. The scheme in question imple- 
ments normalization conditions different from (ITTl) for the theory defined by the full action 
(I7j). No epsilon expansion is employed in this approach, and the only IR fixed point of the 
coupling constant found in four dimensions is stable and trivial, while the three-dimensional 
IR-stable fixed point is nontrivial. 

In summary, our RG analysis recovers, in dimensions D > 2, the decoupling solution 
and one of the scaling solutions [given in Eq. (jlj)] of the DS equations. Moreover, unlike 
the DS equations, the RG analysis is able to show that the decoupling solution is IR-stable 
and hence physically relevant, while the scaling solution (jlj) is IR-unstable. In dimension 
D = 2, the beta function (12T|) is negative for ^ 0, hence the only fixed point, the trivial 
one corresponding to the decoupling solution, becomes IR-unstable. This is consistent with 
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lattice calculations which do not find decoupling (but rather scaling) behavior in the IR at 



D 



111 



What about the other scaling solution of the DS equations given in Eq. ^1 As we 
shall now demonstrate, it can also be obtained from the RG analysis if the so-called horizon 
condition of the original (not refined) Gribov-Zwanziger scenario, 

F{p^) ^ oo for p2 ^ (26) 

in the notation of Eq. ([1]), is implemented as a normalization condition. The horizon con- 
dition is fulfilled for the scaling solutions ^ and ([5]), but not for the decoupling solution. 
Note that we have obtained the scaling solution (jl]) from the RG analysis without ever using 
the horizon condition explicitly. 

Now, in order to implement Eq. fl2Bl) classically, we change the free part ^ of the IR-form 
of the action to 

5{k = / Al m'Al + i d,n-d')d,c'^ + iB'^d.Al^ , (27) 

where the parameter h with the dimension of mass is introduced in order to maintain the 
usual canonical dimensions of the ghost fields. The form of the ghost term in Eq. ( 127|) is 
known as a (isotropic) Lifshitz point in statistical physics 19l. |20|. 

It is clear from Eq. (1271) that the scaling dimensions of the ghost fields differ from the 
usual canonical dimensions, namely, the invariance of the free action under the rescaling 
X — )■ x/s implies the field rescalings 

As for the coupling constants, the three- and four-gluon couplings scale as before in Eq. 
(fT2|) . since the scaling of the gluon field is unchanged with respect to Eq. ( !T0|) . while for the 
ghost-gluon coupling we find 

^?c-Ac^s3-(^/2)(^,^,. (29) 

As a consequence, the upper critical dimension of the theory is now D = Q. We rewrite the 
sole interaction term flT3|) as 

^IR = ~9r'' j d'^x i d.eAl , (30) 
where we have defined a new coupling constant g as g = {b'^/m) gcAc- 
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Following in strict analogy the procedure applied before to the action given by Eqs. ([9]) 
and f|T3l) (details will be provided elsewhere), we again find two fixed points, only that 
this time the trivial fixed point is IR-unstable and the nontrivial fixed point is IR-stable 
in dimensions D < Q. Plugging the value for the stable fixed point into the anomalous 
dimensions yields a scaling solution with IR exponents [in the notation of Eq. ([2])] 

MD) = , «Gp) = ^ . (31) 

Comparing the result to (the approximation to) the solution ^ of the DS equations 
obtained in Refs. , we find that both coincide (exactly) at D = 6, and the difference 
between them in the range 2 < D < 4 is largest for D = 2, where ap = 1/5 (exactly) from 
the solution of the DS equation and ap = 1/6 from Eq. (!3T|) . However, considering the 
usually slow convergence (after Borel resummation) of the epsilon expansion in applications 
to statistical physics, the coincidence of Eq. (ISTi) with the solution of the DS equations 
is surprisingly good. 

Together with the previous dicussion, we now have found two different IR-stable fixed 
points. Which one is physically realized? The realization of the fixed point corresponding to 
the IR exponents fl3Tl) depends on the fine-tuning of the ghost term appearing in the action 
([H]), in order to fulfill the horizon condition. The latter term is perturbatively relevant with 
the scaling dimensions (1251) of the ghost fields. Hence, unless there is a physical reason for 
the fine-tuning of this term or the imposition of the horizon condition (as there actually has 
formerly been thought to be), the fixed point is unstable with respect to perturbations of 
the term in question. 

In summary, the RG analysis also reproduces the second scaling solution ([5]) upon imple- 
menting the horizon condition. However, this condition is IR-unstable under the RG fiow, 
hence in the end the only IR-stable fixed point is the one corresponding to the decoupling 
solution (in dimensions D > 2), in agreement with recent lattice simulations. 
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